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Abstract 

The rate of phase-space exploration in the decay of isovector and isoscalar giant 
quadrupole resonances in 40 Ca is analyzed. The study is based on the time dependence 



S^ of the survival probability and of the spectrum of generalized entropies evaluated in the 

space of lp-lh and 2p-2h states. If the 2p-2h background shows the characteristics typical 
for chaotic systems, the isovector excitation evolves almost statistically while the isoscalar 
excitation remains largely localized, even though it penetrates the whole available phase 
space. 



PACS numbers: 05.45. +b, 24.30.Cz, 24.60.Lz 



Collective nuclear excitations occur on dynamical time scales which are short compared 
to those of the compound nucleus and, therefore, probe simple configurations. In response 
to an external perturbation a collective mode is initially formed as a non-stationary state 
which occupies a small fraction of the available phase space. The subsequent decay, on 
the other hand, involves much longer time scales and explores more complex nuclear 
configurations. The specific characteristics of this process may, of course, depend on 
the energy distribution of the initial state, its multipolarity, its isospin, character or the 
degree of collectivity. The generic features of the time-dependent phase-space exploration 
through decay, however, are of much more general nature |lj] and relate to the quantum 
manifestation of classical chaos in the time-dependent picture |2|, |3], |], [|. The nucleus 
is especially well suited for addressing such questions because of the inherent quantum 
nature and a generic chaoticity of the dynamics. There is also a wealth of experimental 
data which could be useful in verifying some of the theoretical concepts. 

The study of chaos in nuclear physics has been mostly based - so far - on level statistics 
. In practical terms this is rather restrictive since a reliable statistical analysis requires 
very precise energy resolution. It also does not provide firm means for investigating the 
role of collectivity and mechanisms of its coexistence with chaos @]. In this respect the 
study of temporal correlations between an initially prepared non-stationary state and a 
state to which it evolves seems to be much more appropriate. Nuclear giant resonances 
are of central interest in this connection because they are located in a region of high level 
density which is expected to be dominated by chaotic dynamics. To make the theoretical 
studies meaningful one needs a scheme which incorporates the relevant elements, such as 
the possibility of defining a physical collective state, a realistic modeling of the background 
states whose complexity is consistent with the Gaussian orthogonal ensemble (GOE) of 
random matrices |J and, finally, the realistic coupling between the two. 

The recently developed model [HI , based on a diagonalization of the full nuclear Hamil- 



tonian consisting of a mean field part and a residual interaction 

H = Yl € i a l a i + 9 ^ V ij,kl a i a j a l a k, (!) 

i ij,kl 

in the basis of lp-lh and 2p-2h states 

|1) = aja h |0); |2) = o^o^a^o^jO) (2) 

fullfills these requirements and proves numerically manageable ||. A prediagonalization 
of the two-body interaction v in the lp-lh and 2p-2h subspaces defines |1) = XaCf|l) 
and |2) = ^2 C'f 1 2) , and the coupling is mediated by the off-diagonal elements (l|w|2) and 
their complex conjugate. An initially excited state, in response to an external one-body 
field F = J2ij Fij a \ a j can De represented as 

\f) = f\u) = y j f 1 \i). (3) 

1 
The amplitudes Fj = (1|F|0) contain the entire information about the strength of a 
given spectral line corresponding to the state |1) and about the phase coherence among 
these states. On the other hand, the level fluctuations of the states |2) can be used as a 
measure of the degree of complexity of the background states ||. As soon as the coupling 
between the subspaces |1) and |2) is taken into account, the state originally localized in 
the lp-lh subspace, as defined by eq. (^), starts leaking into the 2p-2h space. The degree 
of mixing depends not only on the magnitude of v but also on the nature of the energy 
fluctuations in the 2p-2h space, which can significantly influence the distribution of the 
coupling matrix elements 0, [10| . The most natural quantity for describing the leakage is 



the survival probability, defined as 

P(t) = \(F(0)\F(t)}\\ (4) 

The time evolution of the state vector \F(t)) can be expanded as 

\F(t)) = J2a n e- lEnt/n \n), (5) 



where E n and \n) are the eigenenergies and eigenstates of H in the space of lp-lh and 
2p-2h states. The expansion coefficients are determined by the initial state as a n = 
{n\F(0)). The physical significance of P(t) can be identified from its relation to the 
spectral autocorrelation function G(E) 0, O, |T^| 

P(t) = I dEe- iEt/h G F (E), (6) 

where 

G F (E) = f dE'S F {E')S F {E' + E) (7) 

and S F (E) is the transition strength distribution S F (E) = £ n |F n | 2 <5(.E — E n ). Thus P(t) 
can be obtained from experiments which measure S F (E). Because of finite energy reso- 
lution the experiment determines only an envelope of S F (E) and consequently smoothes 
out the fluctuations in P(t). Average structures which are more interesting are preserved, 
however, provided that the resolution is not too coarse. 

The calculations presented below for quadrupole excitations in A0 Ca are performed 
in the same basis as in ref. ||, i.e. including all lp-lh and 2p-2h states up to 50 MeV 
and using the same residual interaction. We then have 26 lp-lh and 3014 2p-2h states, 
which ensures a realistic description of the transition strength distribution. Since the 
present study concentrates on the phase-space exploration and the role played by chaotic 
dynamics, we distinguish three cases corresponding to different classes of the spectral 
fluctuations in the 2p-2h subspace. As established in ref. || one finds: (a) with no 
residual interaction in the 2p-2h subspace there are many degeneracies in |2) (= |2)) and 
the nearest-neighbor spacing is strongly peaked near zero; (b) inclusion of particle-particle 
and hole-hole matrix elements in (2\v |2') removes all degeneracies and leads to a Poissonian 
distribution of the nearest-neighbor spacings, characteristic of generic integrable systems; 
and (c) use of the full residual interaction yields GOE fluctuations ||, characteristic of 
chaotic dynamics. 

The initial state \F(0)) = F\0) (eq. fl3|)) is already non-stationary in the lp-lh subspace 
and therefore oscillates within the limits set by this subspace. The time evolution of the 



resulting survival probabilities (\F) is normalized to unity) are shown in the upper panels 
of Figs. 1 and 2 for the isovector and isoscalar quadrupole excitations, respectively. It 
is interesting to note that the isoscalar excitation, being more collective, overlaps - on 
average - more frequently with its initial state, as a comparison of the horizontal lines in 
Figs. 1 and 2 indicates. Including the mixing with 2p-2h states (making use of eq. (Ejp) the 
results (Figs. 1 and 2) show that the isovector excitation mixes much more efficiently with 
the background 2p-2h states, both concerning the oscillatory behavior and the asymptotic 
value of the survival probability P(t). The latter systematically decreases with increasing 
degree of complexity in the background states (going from (a) to (c)). Most interestingly, 
for the isovector excitation in the chaotic case (c), P(t) ■ N where N denotes the total 
number (3040) of states in our space, reaches - on average - a value close to 3 (3.08). It is 
known |l|, that a state evolving from generic initial conditions does not visit all the regions 
of the space with equal probability but overlaps more frequently with its initial value. 
This effect, characteristic of quantum ergodicity |lj], is present even if the whole space is 
accessible. In this extreme case it is just the factor of 3 which prescribes the lowest limit 
on the average asymptotic behavior of P(t). Such an 'elastic enhancement' finds empirical 
evidence in nuclear physics |13[ and the factor of 3 is considered as a quantum-mechanical 
signature of chaos [|J. While the factor of 3 is consistent with random- matrix-theory 
estimates |14| the same asymptotic value may, in principle, be associated with the regular 



dynamics if the subspace is defined such that there are no other conserved quantum 
numbers than those defining it. An additional requirement for chaoticity is the initial 
dephasing of P(t) below its asymptotic value |fj, §, ^ |i~2"f . Such a dephasing, indeed, 
takes place for the isovector case, as can be seen from Fig. 3 where P(t) = J dt'P(t')/t 
for all three types of the background spectra is displayed. As a consequence of the high 
density of states the corresponding 'correlation hole' extends over a time interval 
four orders of magnitude longer than the charactistic 'excitation time' of ~ 10 _22 sec. It 
should be emphasized that, contrary to all the previous studies, no ensemble average over 



Hamiltonians or initial conditions is present in our investigations. The isoscalar excitation 
(right panel of Fig. 3) only shows a trace of such a behavior and the asymptotic values 
of P(t) are systematically larger even though the initial state is coupled to the same 
background. 

The question what makes the isoscalar state so strongly localized then arises. Is it a 
manifestation of collectivity or, perhaps, is it that the specific properties of the coupling 
matrix elements block certain regions of the phase space and ergodization only occurs 
in the unblocked regions? A quantity which appears helpful in resolving this question 
originates from the concept of entropy. The information entropy of the state \F) can be 
defined, in a given basis \k), as 

K = -^p fc lnpfc, (8) 

k 

where pk = |(A;|F)| 2 . It provides a quantitative measure of the complexity of the state \F) 
and its localization length in the basis \k) [jn|. The so-defined K is, in principle, basis 
dependent, but the physically preferred basis is determined by the mean field ||16|| . The 
mean field, as the smoothest component of the nuclear Hamiltonian [[HJ, provides a natural 
reference for quantifying the local GOE-type fluctuations. In our case the mean field basis 
corresponds to the unperturbed basis of states |1) and |2). Calculating K{t) along the 
'trajectory' \F(t)) for the isovector and isoscalar states, we obtain asymptotically values 
of 7.10 and 6.53, while the corresponding initial values are 2.56 and 2.40, respectively. 
This is to be compared to K GOE = 7.29 (K GOE = ip(N/2 + l) —0(3/2) [|T§, where ip is the 
digamma function and N is the number of basis vectors). A comparison of these numbers 
indicates non-uniformities in the pk distribution, especially for the isoscalar excitations. 
Actually, even the GOE-type fluctuations result in a gaussian distribution which is non- 
uniform (a uniform distribution maximizes the entropy and corresponds to pk = 1/N 
which, for iV=3040, gives K « 8.02). 

In view of the above mentioned non-uniform phase-space exploration, we find it in- 
structive to calculate the spectrum of g-moments for {pk} and to introduce a generalized 
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entropy |Tj| 

K q = ^~ In f>l. (9) 

1 « k=\ 

From this definition it follows that K qi < K q2 if g 2 < <7i (provided T.kPk = 1)- Equality 
holds for the uniform distribution. For q — > 1 eq. (||) yields the information entropy 
(eq. (||)). The most important property of K q is that with increasing q a higher weight 
is given to the largest components in the set {pk}- For q — ► 0, on the other hand, .fQ 
just counts the number of sites (here the basis vectors \k)) visited, irrespective of how 
frequently they are sampled. For this reason eq. (||) also constitutes a basis for defining 
the multifractal dimensions of non- uniform fractal sets |Tj| . 

For selected q- values Fig. 4 compares the time evolution of K q (t) for the isovector and 
isoscalar excitations when the background states have GOE fluctuations (case (c)). As 
one can see from the large-g behavior of K q (t), which are systematically smaller for the 
isoscalar excitation, the large components of these remain much more localized (larger) 
than those of the isovector excitation. Since, by probability conservation, the number of 
significant components is smaller in the former case, the amplitude of oscillations is larger 
in the corresponding K q (t). On the other hand, the dynamics start to look similar in both 
cases as q decreases and, for q — > 0, K q approaches a value of 8. This signals that, on 
the level of small probabilities, the whole space spanned by 3040 states is visited. This 
aspect of the dynamics is consistent with the scaling properties of the transition strength 
distribution for the isovector and isoscalar states discussed in ref. 0. On the level of 
small components they both scale. The similar calculation of K q (t) for the cases (a) and 
(b) shows that the states evolve to configurations characterized by significantly smaller 
entropies. 

It would be very interesting to verify some of the above predictions experimentally 
- especially the appearance of the correlation hole in the average survival probability 
P{t). This quantity is accessible through the convolution formula (eq. |7D in conjunction 
with eq. (|]). In addition, for the purpose of addressing many specific questions concerning 



quantum-mechanical phase space exploration, the full set of generalized entropies appears 
to be very useful. 
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Figure Captions 

Fig. 1 The time dependence of the isovector quadrupole survival probability P(t) in 40 Ca: 
(lp-lh) no coupling to the 2p-2h subspace; (a) no residual interaction in 2p-2h 
subspace; (b) including only particle-particle and hole-hole matrixelements in the 
diagonalization of the 2p-2h subspace; (c) diagonalization of the full residual inter- 
action in the 2p-2h subspace. The solid horizontal lines indicate the time-averaged 
asymptotic values for the corresponding P(t). 

Fig. 2 Same as Fig. 1 but for the isoscalar quadrupole survival probability. 

Fig. 3 The time-averaged survival probability P(t) in units of 1/N (see text) for the isovec- 
tor (l.h.s.) and isoscalar (r.h.s.) quadrupole resonances: (a) no residual interaction 
in 2p-2h subspace; (b) including only particle-particle and hole-hole matrix elements 
in the diagonalization of the 2p-2h subspace; (c) diagonalization of the full residual 
interaction in the 2p-2h subspace. 

Fig. 4 The time evolution of the generalized entropies defined by eq. (^) for the isovector 
(l.h.s) and isoscalar (r.h.s) giant resonances corresponding to case (c) in Fig. 2. The 
horizontal marks on the left hand side of each panel denote the asymptotic values 
of K q for q = 4 and q = 16. 
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